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ABSTRACT 

We present new near-infrared and optical spectroscopic observations which confirm 
the redshift of the z = 1.44 extremely red object ERO J164502+4626.4 (object # 10 
of Hu & Ridgway 1994; formerly known as 'HR 10' or '[HR94] 10') and a Hubble Space 
Telescope image which reveals a reflected-S~shaped morphology at (rest-frame) near- 
ultraviolet wavelengths. The contrast between the rest-frame far-red (AA8200 — 9800A) 
and near-UV (AA2900 — 3900A) morphologies suggests that the central regions of the 
galaxy are heavily obscured by dust and that the galaxy is most likely an interacting or 
disturbed system. We also present new photometry of this object at 450^m, 850^m and 
1350/um obtained using the SCUBA submillimeter camera on the James Clerk Maxwell 
Telescope. Our sub-mm data are extremely sensitive to emission from cold dust at 
high redshift. The rest-frame spectral energy distribution of ERO J164502-I-4626.4 is 
best understood in terms of a highly reddened stellar population with ongoing star 
formation, as originally suggested by Graham & Dey (1996). The new submillimeter 
data presented here indicate that the remarkable similarity to ultraluminous infrared 
galaxies (ULIRGs) such as Arp 220 and Mrk 231 extends into the rest-frame far-infrared 
which bears the signature of thermal emission from dust, presumably heated by young 
stars. ERO J164502-F4626.4 is extremely luminous {L ^ 7 x lO^'^h^^ ^&) ^^"^ ^^^ty 
(M(just ~ 7 X 10^(T(just/40K)~^/i^Q Mq). If its luminosity is powered by young hot 
stars, then ERO J164502-I-4626.4 is forming stars at the prodigious rate of M = 1000 — 
2000 K^^ Mq yr-^. We conclude that ERO J164502-h4626.4 is a distant analogue of the 
nearby ULIRG population, the more distant or less luminous counterparts of which may 
be missed by even the deepest existing optical surveys. The sub-mm emitters recently 
discovered by deep SCUBA surveys may be galaxies similar to ERO J164502-I-4626.4 
(but perhaps more distant). This population of extremely dusty galaxies may also 
contribute significantly to the cosmic sub-mm background emission. 



Subject headings: cosmology: observations — cosmology: early universe — galax- 
ies: evolution — galaxies: formation — galaxies: starburst — galaxies: individual: 
ERO J164502+4646.4 (HR 10 or [HR94] 10) 



1. Introduction 

Near-infrared (near-IR) imaging surveys have resulted in the discovery of a population of 
infrared-bright, extremely red objects ('EROs'), which may be of significance to studies of galaxy 
evolution. EROs, which we define in this paper as having observed optical-near-IR colors R—K > 6, 
have been identified both in the field, and around high-redshift radio galaxies and quasars (Elston, 
Rieke, & Rieke 1988, 1989, 1991; McCarthy, Persson, & West 1992; Eisenhardt & Dickinson 1992; 
Graham et al. 1994; Hu & Ridgway 1994; Dey, Spinrad, & Dickinson 1995). Images of EROs from 
ground-based telescopes show that they are spatially extended on scales of r-^ O'.'b. This suggests 
that EROs are probably galaxies rather than stellar objects (e.g., brown dwarfs). Most EROs 
were not a priori selected as radio sources and are generally radio quiet. Extremely faint optically 
{R > 24.5), their distances and spectral properties (i.e., whether they are normal or active galaxies) 
remain unknown. 

The extreme colors of these objects might be attributed either to an old stellar population, or 
to a younger stellar population reddened by dust (Graham &: Dey 1996, hereafter GD96; Yamada 
et al. 1997). The spectral energy distributions (hereafter SEDs) of the reddest EROs known are 
inconsistent with unreddened old populations at any redshift, and suggest that they must be highly 
reddened starburst galaxies, perhaps the distant counterparts of the local ultraluminous infrared 
galaxies (ULIRGs; e.g., Soifer et al. 1984) discovered by the Infrared Astronomical Satellite {IRAS 

)■ 

Most EROs have been discovered serendipitously, and their space densities are therefore uncer- 
tain. Some measurements suggest that the most extreme EROs (i.e., those as red as ERO J164502+4626.4) 
are as abundant as quasars, with a surface density in blank fields ~ 0.01 arcmin^^ (Hu & Ridgway 
1994; Cowie et al. 1994). More recent measurements suggest that the surface density of objects 
with R — K' > 6 and K' < 17.5 is ~ O.larcmin^^ (Beckwith et al. 1998) with fainter EROs pos- 
sibly having even higher surface densities (~ 0.7arcmin^^ for EROs with K' < 20; Eisenhardt et 
al. 1998). There is also some evidence that their surface density is higher in regions around high- 
redshift radio galaxies and quasars compared with the general field (e.g., Ar agon- Salamanca et al. 
1994; Dey, Spinrad, &: Dickinson 1995; Yamada et al. 1997). An obvious interpretation of this result 
is that the EROs are clustered around and hence physically associated with the distant, luminous 
AGN, which would imply that EROs are at z ~ 1 — 3 and their rest-frame optical luminosities are 
5-20 times more luminous than unevolved L* ellipticals at the same redshifts. Alternatively, this 
excess number density may suggest that our samples of distant luminous AGN are biased due to 
gravitational lensing by foreground mass concentrations with which the EROs are associated. 

Given that the EROs may constitute a significant population in our Universe, it is important 
to understand their nature. In this paper we present new optical, near-IR and sub-mm observations 
(§ 2) of ERO J164502-F4626.4, an ERO discovered in the field of the z = 3.8 quasar PC1643-F4631 



(object ^ 10 of Hu & Ridgway 1994)^. Our new observations convincingly demonstrate that 
ERO J164502+4626.4 is a z = 1.44 luminous, dusty starburst galaxy, a distant counterpart of 
local ULIRGs like Arp 220 and Mrk 231, as originally suggested by GD96 (§ 3). The detection 
of ERO J164502+4626.4 at sub-mm wavelengths, and the relatively large surface density of the 
ERO population suggests that these objects may be candidates for the sub-mm emitters recently 
discovered in deep SCUBA surveys (Small et al. 1997, Blain et al. 1999a, Hughes et al. 1998, Barger 
et al. 1998), and we briefly discuss this possibility and its consequences in § 3.5. 

We assume Hq = 5O/150 kms^^ Mpc^^, qq = 0.5, and A = throughout. For this cosmology, 
the luminosity distance of ERO J164502-I-4626.4 is 10.53/i^q Gpc. For qQ=0.1, the luminosity 
distance is larger by a factor of 1.286. 



2. Observations and Results 

2.1. HST WFPC2 Imaging 

The field of PC 1643+4631 was observed using the Wide Field Planetary Camera 2 (WFPC2; 
Trauger et al. 1994, Holtzman et al. 1995) on the refurbished Hubble Space Telescope (HST) through 
the F814W filter on UT 1997 May 7. A position angle (PA.) of 68.2° was used in order to 
simultaneously observe both ERO J164502+4626.4 and ERO J164457+4626.0 (objects # 10 and 
14 of Hu & Ridgway 1994). Both objects were targetted on the WF CCDs in order to optimize the 
detection of faint diffuse emission. We obtained 4 exposures over 2 orbits with the telescope dithered 
by ~ 0''7 between orbits. The total integration time was 5300 sec. The images resulting from the 
calibration pipeline were corrected for cosmic rays, registered and coadded. Flux calibration was 
performed using the photometric zero points of Holtzman et al. (1995). The final image reaches a 
3cr limiting surface brightness of 26.77 AB mag arcsec"^ in an in" aperture. 

The WFPC2 F814W image of ERO J164502+4626.4 (rotated to a normal orientation where 
north is up and east is to the left) is shown in figure |l|. The galaxy is extended by roughly 0''9 in 
P. A. ~ —21° in a reflected 'S-shaped' morphology that is suggestive of an interacting system or 
spiral / tidal arms. We spatially coregistered the HST image with the near-IR X-band image of 
GD96 using 16 common objects; the relative positional accuracy is better than 0''02. A comparison 
of the F814W and X-band images reveals that the bulk of the near-IR emission arises from a region 
of low optical flux, in between the two bright lobes of optical emission (the ends of the reflected 
'S'; figure |2|). The brightest region of optical emission is not cospatial with the peak of the near-IR 
emission, but instead lies ~ 0''4 south of it (figure g). 

In a 3" diameter aperture, ERO J164502+4626.4 is found to have a magnitude of 24.6 ±0.1 AB 



^This ERO has been previously referred to as 'HR 10' and '[HR94]' 10 (e.g., Graham & Dey 1996; Cimatti et 
al. 1998). The nomenclature used in this paper follows the lAU recommendations for source names. 



mag (at Aobs ~ 7930A) or, equivalently, a Cousins /-band magnitude of Ic ~ 24.2. The color of 
ERO J164502+4626.4 in this aperture is then Ic — Kk, 5.8=b0.1, consistent with the values quoted 
by GD96 and HR94. The magnitudes of the northern and southern optical lobes, measured in 
apertures of diameter 0''5, are 26.4 it 0.1 AB mag and 26.0 it 0.1 AB mag respectively. 



2.2. Optical Spectrosopy 

We obtained optical spectra of ERO J164502+4626.4 on U.T. 1998 May 4 using the Low 
Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) at the Cassegrain focus of the Keck II 
Telescope. The observations were made in good seeing (FWHM ~ 0"6 — 0''8) through a 1''0 wide 
slit oriented in P. A. = 66.2°. LRIS was configured with the 400 line mm~^ grating (Ab^a^e ~ 8500A; 
AAfwhm ~ lOA), and the spectra cover the wavelength range AA5910 — 9730A. The total exposure 
time was 3 hours, and the target was dithered along the slit after every 30 min sub-exposure in order 
to allow for proper fringe removal and sky subtraction. Wavelength calibration was performed using 
NeAr lamps obtained immediately following the ERO J164502+4626.4 observations, and relative 
flux calibration was performed using observations of the standard stars Wolf 1346 and Feige 34 
(Massey et al. 1988, Massey & Gronwall 1990). Although the observations were made under non- 
photometric conditions, the spectrum was scaled to an absolute scale using published photometry 
(GD96, Hu & Ridgway 1996). 

The sky-subtracted two-dimensional spectrum is shown in figure ^, and the extracted one- 
dimensional spectrum (extracted in a I'll aperture) and the corresponding la errors are presented 
in figure ^ The notable features in the optical spectrum are a strong emission line at A9090A 
(which corresponds to the [OII]AA3726,3729 doublet at z = 1.439) and a very red continuum 
emission {Fy oc iy~5-5±o.5 j^^ ^.j^g wavelength range AA6500 — 9000A). The signal-to- noise ratio in the 
continuum is too poor to convincingly detect any absorption lines or spectral breaks. The measured 
flux, FWHM and rest frame equivalent width of [OII]AA3726,3729 are presented in Table ||. 



2.3. Near-Infrared Spectroscopy 

Near-IR spectroscopy of ERO J164502+4626.4 was obtained using the Cryogenic Spectrograph 
(CRSP; Joyce 1995) at the Cassegrain focus of the KPNO 4m Mayall Telescope on the night of 
U.T. 1997 June 21. CRSP was configured with grating # 4 (200 line mm~^. Ablaze = 3.0 ;um, 
dispersion ~ 12A pixel"^) used in 2nd order and a l''0 slit to provide a resolution FWHM of 0.0042 
/im over the wavelength range 1.444/im to 1.750/im. Observations of ERO J164502+4626.4 were 
obtained by dithering the object along the slit; each dither sequence consisted of five equally spaced 
target positions along the slit separated by 8". The slit was oriented at P.A.=59.3° in order to 
include a nearby galaxy (object 6 in figure |l|; our notation follows that of HR94) to aid in accurately 
coadding the data. Each individual observation was 3 min long, and the total exposure time was 



195 mill. The seeing during the observations was typicahy 0''8 in the K-hand. The atmospheric 
absorption was corrected using observations of the telluric standard HR6064 (GIV) obtained at 
various airmasses both preceding and following the ERO J164502+4626.4 observations. 

The spectral frames were dark-subtracted, flat-fielded, rectified, sky-subtracted, and coadded 
by using object 6 (see figure |^) to determine the accurate relative offsets. The sky subtraction of each 
spectral frame was performed using a 'local' sky constructed from the median of the four preceding 
and four following frames. The summed two-dimensional spectrum is shown in the upper panel 
of figure y. The extracted spectrum was divided by the mean spectrum of the telluric standards 
to correct for the telluric absorption, and corrected for the stellar features by multiplying by the 
mean spectrum of a GOV star constructed by averaging the spectra of HD109358 and HR5868 
from the near-IR spectral atlas of Lancon and Rocca-Volmerange (1992). Since the continuum 
emission from ERO J164502+4626.4 is only marginally detected in these observations, the telluric 
corrections are largely unnecessary, but also provide us with a rough relative flux calibration. 
The resulting spectrum was then scaled by a constant to be consistent with the i:f-band flux of 
ERO J164502+4626.4. 

The calibrated near-IR spectrum of ERO J164502+4626.4 is shown in figure ||, and clearly 
shows the emission line detected by GD96. The line is narrow, and appears to be resolved into 
at least two components. The ratio of the wavelengths of the two brighter components roughly 
matches the ratio of [N II]A6584/Ha, and is consistent with the original identification (by GD96) of 
the line as Ha+[N II]. The extracted spectrum also shows a weak feature at the expected location 
of the [SII]AA6717,6731 doublet. However, given the signal-to-noise ratio of the present data, we 
caution that the [NIIJA6584 and the [SII]AA6717,6731 detections are marginal, and should be treated 
accordingly (e.g., figure ^). 

The emission line measurements are presented in Table |l[ The numbers in this table were 
derived by using the SPECFIT software (Kriss 1994) as implemented in IRAF to fit Gaussians to 
the emission lines in the calibrated unsmoothed spectrum. The Ha and [NH]A6584 emission lines 
were fit jointly by requiring the two features to have equal width and holding the ratio of their 
central wavelengths constant. The new high-resolution near-IR observation presented here provides 
a more accurate measurement of the equivalent widths and FWHMs of the Ha and [Nil] emission 
lines. In these new higher spectral resolution data, the emission lines are found to be narrower and 
have smaller equivalent widths than those derived by GD96. 



2.4. Sub-mm Photometry 

Data were obtained during 1997 June, July and December, and 1998 January and February 
using the Sub-mm Common-User Bolometer Array (SCUBA; Holland et al. 1998). SCUBA has 
two arrays of bolometric detectors which are operated at 0.1 K to achieve sky background-limited 
performance on the telescope at 350-450 and 750-850/im. Three extra bolometers — 'photometric 



pixels' optimized for use at 1100, 1350 and 2000/xm — are positioned around the long- wave array. 

Photometry of sources significantly smaller than the beam are generally performed using the 
central pixels of each array, which are aligned to within an arcsecond of each other, or the pho- 
tometric pixels. The best photometric accuracy is achieved by averaging the source signal over a 
slightly larger area than the beam, so the secondary mirror was 'jiggled' in a filled-square, 9-point 
pattern, covering 4" x 4". During the 'jiggle' the secondary mirror was chopped azimuthally by 
60" at 6.944 Hz. After the first 9 sec 'jiggle', the telescope was nodded azimuthally to the reference 
position (subsequently every 18 sec). At 450^m and 850/xm, we spent 260 min on source, with a 
further 165 min at 1350^m. 

Skydips were performed before, during and after the target measurements to determine the 
zenith opacities, and the telescope pointing accuracy was checked regularly using 1633+382 and 
3C 345. All data were calibrated against observations of Mars and Uranus. The observing conditions 
during the SCUBA runs were excellent, especially during the early 1998 observing periods. 

Reduction of the 450//m and 850/im data consisted of taking the measurements from the central 
bolometer, rejecting spikes, and averaging over 18 sec time intervals. The signal detected by all the 
bolometers is dominated by spatially correlated sky emission (Ivison et al. 1998a), so data from 
the adjacent rings of bolometers were treated in a similar manner. The residual sky background 
was removed using the median of the inner two rings of pixels, after rejecting those with excessive 
noise. This reduced the noise-equivalent flux density at 850/xm to around 90 mJy Hz"^'^, suggesting 
that the effects of rapid sky variability have been removed entirely. Reduction of the 1350/xm data 
followed a similar process, though without the removal of sky noise. 

At 850/xm, we find a clear detection of ERO J164502+4626.4 at a level well above the expected 
extrapolation of the weak radio emission; at 450 and 1350;um, the detection significance is lower, 
but the measured fiux densities are > 3a detections, and are entirely consistent with thermal 
emission from optically thin dust. The fiux density measurements at 450/im, 850;um and 1350/im 
are reported in Table ^. 

The 850/xm and 1350/xm fiux densities reported in this paper are roughly half the values 
reported by Cimatti et al. (1998). Our more extensive dataset was obtained during exceptional ob- 
serving conditions, as evident from the 450/xm detection, and we can exclude a substantially higher 
fiux with some confidence. Furthermore, our data were obtained using a single telescope with 
internally consistent flux calibration, and the 450^m and 850/im photometric data were obtained 
simultaneously. Hence, the relative photometry (and the derived spectral index) are self consistent 
and more accurate than other sub-mm measurements of ERO J164502+4626.4 presented to date. 
There also exist 95/xm and 175/im observations of ERO J164502+4626.4 obtained with the Infrared 
Space Observatory (ISO), which can provide additional consistency checks on the sub-mm obser- 
vations. We will report on these ISO observations in a forthcoming work by Ivison et al. (1999). 



3. Discussion 

3.1. Spectral Properties 

The spectroscopic observations presented in this paper confirm the z = 1.440 redshift of the 
galaxy by resolving the emission line reported by GD96 into two components which match the 
H-alpha + [Nil] identification, the detection of the [Oil] AA3726, 3729 doublet, and the marginal 
detection of the [S11]AA6717,6731 emission lines at the same redshift. 

The Ha/[N11]A6584 ratio (^ 2.4) and the relatively narrow linewidths of the [Oil] and Ha 
emission lines suggest that the ionization is more likely to be due to young, hot stars than an active 
galactic nucleus (AGN). The spectrum of ERO J164502+4626.4 is similar to that observed in local 
star-forming galaxies (e.g., Kennicutt 1992). The FWHM (deconvolved FWHM ^ 260 km s^^) of 
the [Oil] doublet and equivalent width of the Ha emission line from ERO J164502+4626.4 are 
more typical of star-forming galaxies than AGN (e.g., Liu & Kennicutt 1995). However, since 
an AGN could be hidden behind dust and rendered invisible at rest frame optical and ultraviolet 
wavelengths, a firm statement on the nature of the energy source must await future mid- and far- 
infrared spectroscopy. In the following discussion, we assume that young stars power the entire 
far-infrared emission from ERO J164502-I-4626.4. 



3.2. Morphology 

The HST WFPC2 image (figures |l] and |) clearly shows that ERO J164502-h4626.4 has an 
elongated and distorted morphology. At a redshift of 1.44, the WFPC2 F814W filter samples 
the rest-frame wavelength range AA2900— 3900 (i.e., roughly rest-frame [/-band light), which in- 
cludes the [O II]AA3726,3729 emission doublet. The equivalent width of this emission feature in 
ERO J164502-I-4626.4 implies that the emission line contamination in the F814W filter is ^8%; 
the observed morphology is therefore dominated by the continuum emission from the galaxy. Since 
the filter samples the continuum at near-UV wavelengths just shortward of the 4000A break, the 
morphology observed in the HST image only reflects the distribution of UV bright populations 
(i.e., hot, young stars), and in addition may be heavily modified by the distribution of dust in the 
galaxy. The peculiar UV morphology may be due to spiral arms, tidal features or something more 
complex, and suggests that ERO J164502-I-4626.4 may be an interacting or distorted system. 

The most striking aspect of the rest frame near-UV morphology is that it is very different from 
the rest frame far-red morphology observed in the X-band: the red emission is more symmetric (at 
least at the 0''6 resolution of the GD96 ground-based image) and peaks in a region where the near- 
UV emission falls to a minimum, approximately at the center of the reflected S-shaped structure 
seen in the WFPC2 image (figure ^). Since the far-red emission in galaxies is generally dominated 
by late type giant stars and an older main sequence population and the near-UV emission is 
dominated by younger, hotter stars, the observed morphological difference could be partly due to 



spatially distinct populations, with the older stars being more centrally concentrated (perhaps in 
a bulge?) compared to the younger stars (perhaps in a disk or tidal arms?). 

Is it possible that the UV morphology is dominated by young super star clusters? The specific 
luminosities of the northern and southern lobes observed in the near-UV {Xrest ~ 3250A) image 
are 5.4 x 10^^/i^g ergs~^ Hz~^ and 7.9 x 10^^/i^q ergs~^ Hz~^ respectively. Therefore, these regions 
have luminosities that are at least 3 times larger than the brightest of the super star clusters 
observed in local starburst galaxies (e.g., Meurer et al. 1995), and it is possible that these regions 
are composed of groups of star clusters. We have not corrected the observed flux densities for the 
effects of extinction due to dust, and the luminosities quoted should be regarded as lower limits. 

However, since the extreme red color of ERO J164502+4626.4 requires a signiflcant amount of 
dust, the observed morphological difference could also be produced by a spatially non-uniform dust 
distribution which reddens the central regions more than the outer structure. Indeed, the colors in 
any resolution element are redder than any reasonable stellar population, implying that the entire 
visible spatial extent of the galaxy is strongly reddened. It is intriguing that similar UV/optical 
morphological differences are observed in the inner regions of some nearby ULIRGs where they are 
commonly interpreted as resulting from spatial differences in the dust extinction. 



3.3. Dust 

The SED of ERO J164502+4626.4 is remarkably similar to those of low redshift ultraluminous 
galaxies like Arp 220 and Mrk 231 (figure |6|), and this strongly suggests that ERO J164502+4626.4 
is a dust-enshrouded system. The measured sub-millimeter fiux densities at 450/im, 850//m and 
1350^m are greatly in excess of the power-law extrapolation of the non-thermal radio emission, and 
it is therefore safe to assume that the sub-millimeter flux results from thermal emission from dust 
heated by an AGN or in star forming regions. The mm/sub-mm spectral index is also consistent 
with this hypothesis. 

If we assume that the SED of ERO J164502+4626.4 is indeed identical to that of the local 
ULIRGs, we can estimate the luminosity and mass of the warm dust responsible for the thermal 
emission, and attempt to constrain the star-formation rates necessary to power this galaxy. Using 
the template SED of an luminous IRAS galaxy compiled by Guiderdoni et al. (1998), we find that 
a least squares fit of the Guiderdoni SED to the three sub-mm points give a bolometric luminosity 
of L = 7 =b 1 X 10^^ h^Q Lq with reduced x^ = l-l (figure |^. Almost all of this emission emerges in 
the far infrared: LpiR ~ 6.7 x 10^^ /i^q Lq, where we have defined LpiR as the luminosity emerging 
between rest wavelengths of 10/im and 2 cm. This ranks ERO J164502+4626.4 among the most 
luminous (unlensed) IR galaxies known, and classifies it as an ULIRG. 

Note that the above estimate of the bolometric luminosity requires a large extrapolation. The 
sub-mm luminosity (estimated from frestLJ^*^**[850^m/(l + z)], i.e., the best-determined sub-mm 
flux) is 6 X 10^^ Lq, and so this is an extrapolation of two orders of magnitude. However, the 
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sub-mni flux densities are well-fit by optically thin thermal dust emission with an emissivity law 
with index 1.5 (i.e., Qa oc A"-*^'^) and temperature of T^ust = 40 ± 4 K. This dust temperature is 
typical of ULIRGs (Klaas & Elsasser 1993). These modified black body model fits to the sub-mm 
data have steeper spectra than the Guiderdoni et al. SED at Arest ^ 40 //m, and result in slightly 
lower values of Lfir. However, although the sub-mm data are all on the Rayleigh- Jeans tail of 
the dust emisssion and therefore do not strongly constrain the dust temperature, the small x^ 
of the fit and the reasonable value of the derived dust temperature suggest that our estimate of 
the far-infrared luminosity is uncertain by at most a factor of two. The corresponding mass of 
dust is Mfiust ~ 7 X 10^(T(just/40K)~^/i^Q Mq for optically thin emission from grains with normal 
interstellar parameters (Draine & Lee 1984, Hildebrand 1983). 



3.4. Star Formation Rates 

The star formation rates can be estimated from either the Ha luminosity or the sub-mm 
continuum luminosity only if we assume that young, hot stars provide the energy source for the 
line and thermal continuum emission. 

For our adopted cosmology (i7o=50/i5o kms^^ Mpc^^ and go=0.5), the luminosity in the Ha 
emission line is Lho ~ 4.4 x 10^^ /i^q erg s^^. If this emission line is powered entirely by young, 
hot stars, then its luminosity implies a total star-formation rate of ~ 40 /igg Mq yr^"*^ (Kennicutt 
1983). This is almost certainly a lower limit to the true star formation rate in the galaxy, since we 
have not made any corrections for the dust extinction. 

A more reliable estimate for the star formation rate in FRO J164502-I-4626.4 can be derived 
under the assumption that the sub-mm continuum excess is due to thermal emission from dust 
grains heated by young, hot stars. If stars are formed at a constant rate over 10^ — 10® yr with 
a Salpeter initial mass function {4>{m) oc m"^'^^; 0.1 < m < 100 Mq), a luminosity of 10^^ Lq 
corresponds to a star formation rate M between 14 — 24M0yr^"'^ (Leitherer and Heckman 1995). 
Therefore the observed bolometric luminosity corresponds to a total star formation rate of M = 
1000 — 1800 /i^Q Mq yr~^ (or a factor of three smaller for the formation rate of only massive stars). 

A different estimate for the formation rate of massive stars can be derived by relating the 
total luminosity to the mass consumption rate in early type stars. Scoville & Young (1983) derive 
MoBA ~ 7.7x W-^^ (L/ Lq) Mq yr-\ which implies a formation rate of ^ 560(rdust/40K)5 Mq yr'^ 
for the massive stars. Note that this estimate does not account for the formation of lower mass 
stars, the mass cycled through or locked up in lower mass stars and their remnants. This star 
formation rate is therefore a lower limit, and as such is consistent with the estimates derived above 
using the Leitherer & Heckman (1983) starburst models. 

If the sub-mm luminosity is indeed powered by young stars, the derived star formation rate 
is 26 — 45 times greater than that inferred from the Ha luminosity; if this difference is due to 
dust extinction by a foreground screen, it implies an extinction of Ay ~ 4.5 mag. For a Galactic 
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extinction curve, this would imply that the intrinsic ratio of [O II] to Ha of (-f'[oiil/-^HQ!)o ~ 1-3, 
which is at the extreme limit of observed values for local star forming galaxies (e.g., Kennicutt 
1992). Since the dust is more likely to be intermixed with the line-emitting gas than external 
to it, the extinction derived by comparing the sub-mm and Ha luminosities may not be directly 
applicable to the reddening correction of the emission line ratios. The Ay inferred here is about 
2.5 mag larger than that estimated by fitting the ERO J164502+4626.4 optical and near-IR flux 
densities with a reddened Sb galaxy model SED (GD96). 

The star formation rate derived from the sub-mm continuum emission is extremely large and 
may lead us to the conclusion that the dust emission is powered by an AGN rather than by star 
formation. Indeed, the source of the large luminosities in local ULIRGs is still a matter of debate 
(e.g., Sanders & Mirabel 1996). However, there is no evidence from the existing spectroscopy that 
this is the case: as noted in § 3.1, the emission line ratios and widths are more typical of starburst 
galaxies than AGN. 

The dust mass estimate from the previous section and the star formation rate derived here 
can be combined to yield a crude order-of-magnitude estimate for the lifetime of the starburst. 
Assuming that ERO J164502+4626.4 has a 'normal' gas-to-dust ratio of 0.01, the total mass of gas 
is ~ 7 X 10"'^'^ Mq, and the lifetime of the burst is iburst ~ 1-3 x 10^ yr. The total mass of stars 
produced is ~ 10^^ Mq, comparable to that of a present-day massive galaxy. 

In this section we have assumed that the sub-mm continuum emission is due to thermal emission 
from dust heated by young, hot stars. One curious property of ERO J164502+4626.4 which may 
be relevant to this hypothesis is that the rest frame 60//m emission predicted by the dust emission 
and the nonthermal radio emission at a rest wavelength of 6 cm (extrapolated from the flux density 
observed at 1.4 GHz assuming a S^ cc i/"*^'^ synchrotron spectrum) does not follow the well-known 
60//m - 6 cm correlation obeyed by local star forming galaxies (e.g., de Jong et al. 1985). The rest 
frame 6 cm emission predicted by this correlation for ERO J164502+4626.4 is ~ 300/iJy, more than 
an order of magnitude above the detected radio emission (Table |2|) . The only ways of decreasing the 
predicted 60/im emission would be to decrease the dust temperature and increase the dust optical 
depth, both of which result in worse fits to the sub-mm data. The existing ISO data may be able 
to help constrain the dust temperature and optical depth estimates (Ivison et al. 1999). This issue 
can also be resolved with better far-infrared observations of ERO J164502+4626.4, which will soon 
be possible with the Space Infrared Facility (SIRTF). 

If the large deviation from the 60//m - 6 cm correlation is real, its interpretation is unclear. 
It is unlikely that this deviation is due to the sub-mm flux being heated by an AGN instead of a 
starburst, since in the local universe AGN tend to depart from the correlation by having excess 
nonthermal radio emission, not excess far-infrared emission (e.g., Dey &: van Breugel 1994). Some 
local ULIRGs (e.g., Arp 220) do show low frequency turnovers in their radio spectra which are 
generally attributed to thermal absorption of the radio emission (Sopp & Alexander 1991). On the 
other hand, the first reliably identified sub-mm selected high redshift galaxy, SMM J02399— 0136 
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(z = 2.80; Ivison et al. 1998b), has a 'normal' 60/im - 6 cm ratio, yet shows some evidence for an 
AGN contribution to the total emission. Since the far-infrared flux densities are poorly constrained 
at present, and the origin of the correlation in local galaxies is also not completely understood, we 
hesitate to draw any firm conclusions from this observation. 



3.5. Cosmological Implications 

Measurements of the global star-formation history of the Universe, using deep redshift surveys, 
e.g., the Canada France Redshift Survey (Lilly et al. 1996) reaching z ~ 1, the statistics of Lyman- 
limit galaxies (Steidel & Hamilton 1992) at z = 3.4, the Hubble Deep Field (HDF) 2.5 < z < 4 
(Madau et al. 1996, Connolly et al. 1997), imply that the star-formation and metal-production 
rates were about 10 times greater at z ~ 1 than in the local Universe, that they reach a maximum 
somewhere in the redshift range 1 < z < 1.5 and remain roughly constant (or perhaps slowly 
decline) at higher redshifts (Steidel et al. 1999). 

These conclusions, which are based almost entirely on samples selected at optical and near-UV 
wavelengths, may be misleading (Small, Ivison & Blain 1997; Blain et al. 1999a). Absorption by 
dust in regions of star-formation may have distorted our picture of galaxy evolution in the high- 
and low-redshift Universe in two ways. First, neglect of dust leads to an underestimate of the the 
star-formation rate in known high- and low-redshift objects. Second, it is possible that an entire 
population of heavily dust-enshrouded high-redshift objects, such as EROs, has escaped undetected 
in the optical/UV surveys. 

A useful diagnostic of the redness of the SED is the rest frame infrared-to-blue luminosity ratio 
(Lpjji/Lb), where LpiR = vL^ at Arest = 80/xm and Lb ^ I'Liy at Arest = 4400A. By this measure, 
the reddest galaxy in the UGC is Arp 220 {L fir/ Lb ~ 60; Soifer et al. 1984). The mean Lfir/Lb 
value for the ultralumimous sample is 25 (Sanders et al. 1987), and the reddest ultraluminous 
galaxy in the IRAS Bright Galaxy Catalog is IRAS 12112-1-0305, which has Lfir/Lb ~ 70. If 
we estimate Lb for ERO J164502-I-4626.4 by interpolating the observed / and J fluxes, we find 
that ERO J164502-F4626.4 has Jb ~ 2.3 /iJy. Hence, Lb « i^L^(4400A) « 2.2 x IO^^Lq, and 
Lfir/Lb ~ 300! The name ERO is therefore justified for this remarkable object. 

The relevance of ERO J164502-I-4626.4 to cosmology and galaxy formation depends critically 
on the space density of the ERO population, which is highly uncertain at present. Under the as- 
sumptions that the EROs form a homogeneous population and the redshift and sub-mm luminosity 
of ERO J164502-I-4626.4 is typical of this class, GD96 estimated that the space density of these 
objects is pero ~ 2.1 x IO^^/i^q Mpc~^ mag'^. Hence, EROs may be more abundant than the 
IRAS ultraluminous galaxies and quasars by two orders of magnitude. If the star formation rate 
for these galaxies is similar to that observed in ERO J164502-I-4626.4 (i.e., M > 580 Mq yr^^), the 
star formation rate at z = 1.4 associated with the ERO population is ^ OAh^o MQyr~^ Mpc~^. 
This sub-mm estimate of the star formation density exceeds the rest-frame UV estimates of the 
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high-redshift star formation rate (from the CFRS and HDF surveys) by more than a factor of six. 
It is noteworthy that this crude estimate is roughly consistent with the predictions of Pei & Fah 
(1995) (see also Hughes et al. 1998). 

Deep surveys at 850/im of the sub-mm sky have resulted in the detection of several faint sources 
(^ 4 mJy), about 20% of which appear to have no obvious optical counterpart brighter than I ^^ 25 
(Smail et al. 1998). It is possible that these 'missing' optical identifications are galaxies similar 
to ERO J164502+4626.4, perhaps even fainter at optical wavelengths, but with similarly strong 
thermal dust emission at sub-mm wavelengths. Since the cumulative surface density of the sub-mm 
continuum emitters with 850/im flux densities > 4 mJy is 2.4 it 1.0 x 10^ deg~^, the surface density 
of the optically faint fraction is ~ 480 deg~^ (estimates here are based on the surface density 
determined by Blain et al. 1999b; however, see Barger et al. 1998 for a lower estimate). The surface 
density of EROs is highly uncertain: Hu & Ridgway (1994) estimated that the surface density of 
objects with I — K > 6 and X < 19 is ~ 36 deg~^, whereas estimates of the surface density of less 
extreme EROs (with R — K > 6) discovered serendipitously tend to be larger by at least an order of 
magnitude (e.g., Dey, Spinrad & Dickinson 1995; Knopp & Chambers 1997; Beckwith et al. 1998). 
The roughly comparable space densities of these two populations suggests that it is possible that 
EROs comprise a significant fraction, perhaps all, of the optically faint sub-mm emitters discovered 
in recent surveys. A large fraction of both local ULIRGs and distant SCUBA sources (Smail et 
al. 1998) appear to be interacting or merging systems, perhaps similar to ERO J164502+4626.4. 

Since the sub-mm sources thus far detected at flux densities > 0.5 mJy arguably account for 
upto 100% of the 850/im cosmic sub-mm background (Smail et al. 1997, Blain et al. 1999), the 
ERO population may comprise a substantial component of the resolved emission. The relevance of 
the properties of ERO J164502+4626.4 for our general understanding of both the ERO and sub- 
mm populations will hinge on whether it represents a particularly high luminosity or low redshift 
member of this group. 

The very red SED of ERO J164502+4626.4 also implies that it will be extremely difficult 
to detect the most actively star-forming galaxies at higher redshifts, since they may be com- 
pletely enshrouded by dust and undetectable at optical and near-IR wavelengths. For instance, 
if ERO J164502+4626.4 were at z = 5, it would have a ET w 25.8 and / > 32 ( Ho=50, qo=0.5), and 
therefore undetectable in even the deepest existing ground- and space-based surveys. In contrast, 
the large negative sub-mm k-correction would imply that the 850^m flux density of such an object 
would be ~ 3 mJy, comparable to that at z = 1.44. The existence of dusty star- forming systems 
at 2 > 5 is not an unrealistic expectation: the z = 5.34 galaxy 0140+326RD1 (Dey et al. 1998) 
has been recently argued to be a dusty star-forming system based on its red I — J color (Armus et 
al. 1998). It is important to note, however, that the inferred dust content of 0140+326RD1 is much 
lower than that of ERO J164502+4626.4: the implied extinction is Ay ~ 0.5 mag for 0140+326RD1, 
compared with Ay ~ 4.5 for ERO J164502+4626.4. Investigating the properties of the higher red- 
shift counterparts of systems like ERO J164502+4626.4 will require sensitive ground-based sub-mm 
facilities with high-angular resolution and space-based telescopes operating at near-IR and mid-IR 
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wavelengths. In the near-future, the Space Infrared Telescope Facility {SIRTF) will provide flux 
density observations of EROs at mid- and far-infrared wavelengths (~ 3 — 180 jiva) allowing us 
to place more critical constraints on the dust masses, temperatures and source luminosities, and 
perhaps determine the relationship of these objects to present-day galaxies. 



4. Conclusion 

We have presented new optical, near-IR and sub-mm observations of ERO J164502+4626.4 
(object # 10 in Hu & Ridgway 1994), an extremely red object in the field of the QSO PC 1643+4631A. 
The new optical and near-IR spectroscopy show that ERO J164502+4626.4 is a distant galaxy ly- 
ing at a redshift of z = 1.44. The peculiar rest frame near-UV and far-red morphologies suggest 
that ERO J164502+4626.4 is a disturbed or interacting system. The far~IR and sub-mm evidence 
are consistent with the hypothesis that this ERO is a dust-enshrouded object, with its luminosity 
(L Ri 7 X 10^^ h^Q Lq) powered by either a starburst or an AGN. The existing spectral data for 
ERO J164502+4626.4 show no strong evidence for AGN emission at rest-frame near-UV or optical 
wavelengths. If the observed sub-mm continuum flux is due to optically thin thermal emission from 
dust heated by a young, star forming population, the inferred star formation rate is extremely large 
(1000 - 2000 MQyr-'^). In all its known properties, ERO J164502-F4626.4 appears to be a high- 
redshift luminous counterpart of the dusty, ultraluminous galaxies discovered in the local Universe 
by IRAS. Although it is difficult to draw general conclusions from a single object, less luminous 
or higher redshift galaxies similar to ERO J164502-I-4626.4 may be missed from even the deepest 
existing optical and surveys. It is important to determine the space densities of these EROs and 
their relevance to our understanding of galaxy formation and evolution in the distant universe. 
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WFPC2 F814W 




Fig. 1.— HSTWFFC2 image of ERO J164502+4626.4 (object 10 in Hu k Ridgway 1994) obtained 
through the F814W filter. The field of view shown is 50" on a side, and the ERO and the QSO 
PC1643+4631A are labelled. North is up and east is to the left. ERO J164502+4626.4 ("HR 10")is 
located at a = 16^45^02536, 6 = +46°26'25^'5 (J2000), and the offset from the QSO is Aa = 
+13''85, AS = +10''27. The galaxy labelled "6" was used for aligning the dithered IR spectra (see 
text). The inset, 4" on a side, shows ERO J164502+4626.4. 



19 



F814W 






1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



o 




\3> 



"l I I I I I I I I I ll2 L_2_l I I I L 




-2 2 

Aa{") 







Fig. 2. — This figure presents a comparison of the observed K and F814W morphologies of 
ERO J164502+4626.4. The left panel shows a contour plot of the HST WFPC2 image obtained 
through the F814W filter. The contour levels are drawn at levels (2,3,4,5,6,7,8,9,10) xa^fcy, where 
asky = 25.56 ABraag arcsec"^. Negative contours are represented by dotted lines. The right panel 
shows a contour plot of the Lucy-deconvolved X-band image from Graham & Dey (1996). The 
contour levels are drawn as in the left panel, but with asky equivalent to 21.38 mag arcsec^^. The 
resolution of the F814W image is ~ 0''1 whereas the resolution of the Lucy-deconvolved X-band 
image is ~ 0''28. Note that the morphologies are different and that the peak of the X-band emission 
occurs in a region of low optical emission. 
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Fig. 3. — The above panels present portions of the two-dimensional spectra of ERO J164502-I-4626.4 
showing the emission line detections. The upper panel shows the detection of the Ha emission line in 
the coadded near-IR spectrum obtained using CRSP on the KPNO 4-m telescope. The lower panel 
shows the detection of the [Oil] AA3726, 3729 emission doublet in the optical spectrum obtained 
using LRIS on the Keck II Telescope. 
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Fig. 4. — Optical spectrum of ERO J164502+4626.4 obtained using LRIS on the Keck II Telescope. 
The upper panel shows the observed spectrum extracted in a l''0 x l''7 aperture smoothed using 
a 13A width boxcar filter, and the bottom panel shows the corresponding la error spectrum. The 
[OII]AA3726,3729 emission doublet is marked. 
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Fig. 5.— Near-IR spectrum of ERO J164502+4626.4 obtained using CRSP on the KPNO 4- 
m telescope. The upper panel shows the observed spectrum extracted in a l''0 x l''7 aperture 
and smoothed using a 36A width boxcar filter (solid line), along with the relative atmospheric 
transmission (long dashed line). The bottom panel shows the corresponding la error spectrum 
which is dominated by the OH telluric emission lines. The Ha emission line and locations of the 
[NHJA6584 and [Sn]AA6717,6731 emission doublet are marked. 
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Fig. 6. — The SED of ERO J164502+4626.4 between the radio and optical wavebands, repre- 
sented by filled circles. The right-hand scale gives flux densities for ERO J164502-I-4626.4. For 
comparison, we have plotted the SEDs of the ultraluminous IRAS sources F 10214-1-4724 (Rowan- 
Robinson et al. 1993; Barvainis et al. 1995), Mrk 231 and Arp 220 (D.H. Hughes priv. comm.), 
and SMM 02399—0136 (Ivison et al. 1998b) with units of luminosity density (left-hand scale). 
These lines are broken in regions where only upper limits are available. For F 10214-1-4724 and 
SMM 02399—0136 the SEDs are corrected for lensing by factors of 30 and 2.5 respectively (e.g., 
Graham & Liu 1995; Broadhurst & Lehar 1996; Ivison et al. 1998b). 
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Fig. 7. — The rest frame far-infrared SED of ERO J164502+4626.4 wavebands (filled circles) 
compared with a Guiderdoni et al. (1998) semi-empirical SED for IR luminous galaxies and a 
modified blackbody {Fy oc By[l — exp{— (z//fo)'^}]) of temperature T = 40 K and emissivity index 
(3 = 1.5. The models shown were fitted to the three sub-mm flux density measurements. 
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Table 1. Emission Line Measurements of ERO J164502+4626.4^ 



Line 



^obs 



Redshift 



Flux 



FWHM M/A.rcst' 



[OII]AA3726,3729 9090.6±0.9 

Hq 16030±7 

[NII]A6584 16081.1^ 

[SII]AA6717,6731 16437±7 



1.439 


2.6±0.4 


420±100 


47±5 


1.443 


33±8 


597±140 


89±20 




14±5 


597^ 


37±13 


1.444 


14±6 


532±250 


38±17 



^All quoted measurements are based on Gaussian fits to the emission lines. The spectra 
were obtained through l!'0 wide slits, and the spectral extractions used in these measure- 
ments are l"7 wide in P.A.=66.2 ([Oil]) and P.A.=59.3 (Ha, [Nil] and [SII]). 

Rest-frame equivalent widths assume z = 1.440. 

■^The central wavelength and width of the [NII]A6584 emission line are fixed with respect 
to the derived values for the Ha emission line. 



Table 2. Photometry of ERO J164502+4626.4 



Observed 


Rest 


Flux Density 


Detector/ 


Reference 


Wavelength 


Wavelength 




Instrument 




4400A 


I8OOA 


0.16 ± 0.07 /iJy 




GD96, Hu & Ridgway (1992) 


7930A 


3250A 


0.52 ± 0.06 /iJy 


WFPC2/i/5'r 


This paper 


1.2/im 


4920A 


6.4 ±2.1 /liJy 




GD96, Hu & Ridgway (1992) 


1.6/xm 


6560A 


14.8 ± 3.6 /iJy 




GD96, Hu & Ridgway (1992) 


2.2^m 


9OIOA 


27.7 ± 0.6 ^Jy 




GD96 


12/im 


4.9/im 


< 75 mjy 


IRAS 




25/jm 


10.2^m 


< 60 mJy 


IRAS 




60^m 


24.6^m 


< 84 mJy 


IRAS 




lOO^m 


41.0^m 


< 270 mJy 


IRAS 




450^m 


184^m 


32.3 ± 8.5 mJy 


SCUBA/JCMT 


This paper 


850^m 


348^m 


4.89 ± 0.74 mJy 


SCUBA/JCMT 


This paper 


1350//m 


553/im 


2.13 ±0.63 mJy 


SCUBA/JCMT 


This paper 


3.6 cm 


1.5 cm 


35 ± 11 ^Jy 




Frayer (1996) 


20 cm 


8.6 cm 


< 300 ^Jy 




Frayer (1996) 



